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EXTENSION of bodhdahy -layer head? -transfer theory to 
COOLED TURBINE BLADES ' 

By W. Byron Brown and. Patrick L. Donougbe 


SIM4AEY 

A review of existing Bound ary -layer heat -transfer theory appli- 
cable to outside heat -transfer coefficients was made. The influences 
of Mach number , temperature ratio, and exponents of gas -property 
temperature relations were computed and were shown to be relatively 
small, -when the gas is air, for Mach numbers less than 2 with tem- 
perature ratios between 1 and 4. An equation for the average heat 
transfer of a surface was derived for constant wall temperature when 
the boundary layer changed from laminar to turbulent on the surface. 
These results indicated that the parameters needed to calculate the 
Nusselt number are the Reynolds and Euler numbers and the transition 
ratio (length of surface with laminar boundary layer divided by total 
length of surface ) . Good agreement between the average heat transfer 
predicted by the theory and same experimental results from cylinders, 
an airfoil, and turbine blades was obtained. 


UJTROEJUCTI0N 

Mo satisfactory coir elation of heat transfer from turbine blades 
is possible between cooled turbine blades of the impulse and the 
reaction type when the correlation schemes used for boiler pipes, 
heat exchangers, and reciprocating engines are applied. It is known 
that the convective heat-transfer rates between the fluids and the 
metals involved are dependent on the boundary layers. In the general 
case of heat transfer through the boundary layer, the stream tempera- 
ture may differ appreciably frcm the wall or surface temperature, so 
that variations of specific heat, density, viscosity, and thermal 
conductivity should be considered, as well as the effects of Mach 
number and pressure gradient in the direction of flow. 

Theoretical studies have been made of heat tr ans fer through 
laminar and turbulent boundary layers of certain simple types. One 
of the earliest laminar studies was made by Pohlhausen (reference l) , 
who assumed no pressure gradient (flat plate) , various Prandtl num- 
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bers, and constant fluid properties through the "boundary layer. This 
analysis was extended to include large Mach numbers (supersonic flow) 
and temperature ratios (the absolute temperature of the fluid outside 
the boundary layer divided by the absolute temperature of the wall) 
widely different from unity (references 2 to 4 ) . Variations In 
specific heat and the temperature -viscosity law are. also considered 
in reference 4. In. references 5 and 6 , the pressure gradient along 
the flow path is included, but constant fluid properties are assumed, 
thus Implying temperature ratios near 1. All these analyses assume 
a constant wall temperature. Reference 7 considers variations in the 
wall temperature without a pressure gradient in the flow direction. 

In the case of theoretical formulas for turbulent boundary 
layers, Reynolds analogy is usually employed. Summaries of the 
resulting formulas for heat transfer are given In references 8 and 9. 
The problem for incompressible flow In a pipe and over a flat plate 
is discussed in reference 10. 

Same experimental results have been published on average outside 
heat -transfer rates between the gas and the turbine blades from sta- 
tionary cascades. Eo nozzles preceded these cascades. Data on 
impulse blades with temperature ratios of 0.8 and 0.9 are given in 
references 11 and 12, respectively, and reference 13 gives results 
for temperature ratios above 1. Even though these results are 
plotted in the same manner, different correlations are shown for dif- 
ferent cascades (reference 14). In order to unify such results, 
theoretical studies are needed and new parameters not used in tube 
investigations may be necessary. 

The theoretical effects of the various parameters on the heat 
transfer with a laminar boundary layer and constant wall temperature 
were analyzed at the KACA Lewis laboratory and are discussed herein. 
"With the assumptions of constant wall temperature and Mach numbers 
less than 2, a formula is derived for the average outside heat- 
transfer coefficient at points where the boundary layer Is partly 
laminar and partly turbulent. The effect of two additional param- 
eters, the Euler number and the transition ratio, are calculated 
and presented. Results from experimental investigations are assembled 
and compared with results calculated by means of the theoretical 
formulas. 
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STATUS OF BDOHDAHT-IAXEER TEEGRT TOTH HEAT TRAUSEER 

When a fluid flows along a solid 'boundary, -the part in. immediate 
contact with the surface has zero Telocity relative to the solid 
(reference 15, pp. 676-680); whereas in the free stream, the fluid 
has its greatest velocity relative to the solid. Between the wall and 
the free stream, the velocity changes continuously from zero to the 
free -stream value; this region is called the boundary layer. One of 
the characteristic features of this region is that it is very thin 
compared with the length of the "body for the ranges of Reynolds num- 
ber’s considered herein (10^ < Re < 10 6 )-. Under these circumstances 
(thin boundary layer) , the same equations hold for curved surfaces 
as for plane ones for a laminar boundary layer (reference 15, 
p. 120) , so that the radius of curvature will not enter into the 
solution in general (the stagnation point is an exception) . Such a 
boundary layer is shown in figure 1, which greatly exaggerates the 
relative thickness of the boundary layer. The solid surface is shown 
in the figure as a plate placed in the middle of a converging channel 
and, therefore subject to a decreasing pressure gradient for subsonic 
flow. The distance x is measured along the surface from the stag- 
nation point and the y direction is perpendicular to the surface. 
(All symbols used in the report are defined in Appendix A. ) 


Equation for Heat -Transfer Coefficient 

At the wall surface, heat leaves the fluid by conduction at the 
rate per unit time per unit area k^CcT/Sy)^., where is the 

thermal conductivity of the fluid at the wall temperature and. 

(dr/&y) v is the temperature gradient perpendicular to the wall at 

the wall. This rate of conduction is equal to the amount of heat 
that enters the wall per unit time per unit area, usually given as 
the product of a heat -transfer coefficient H and a temperature dif- 
ference between the fluid in the main stream and the wall . When the 
fluid is moving rapidly, an effective temperature T is used in the 
heat -transfer equation. Thus 


H(T e 




( 1 ) 


T is not a function of x. 
w 


where 
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Tile effective temperature T e is not the temperature indicated 

by a thermometer moving with the fluid, nor is it the stagnation tem- 
perature indicated by a thermometer placed in a calming chamber where 
the fluid is brought to rest without heat transfer; that is, adiabat- 
ically. It is instead the temperature that satisfies equation (l) at 
a specific point x when (dT/dy}^. approaches zero and E remains 

constant. Thus, T = T in this case and, consequently, T is 

often called the adiabatic wall temperature. By utilizing this con- 
cept for the effective gas temperature, the differential equations of 
the boundary layer can be solved for (c*T/c>y) and T , and H may- 
be computed, through the use of equation (l) . 


Equations of Laminar Boundary Layer 

The equations of the compressible laminar boundary layer for 
steady-state flow of an actual gas with heat transfer are (refer- 
ence 4); 

Momentum equation, 

(2) 

Continuity equation, 

(pu) + 5^ * 0 ( 3 ) 

Energy equation, 

j °p ( pu M + pT If) ■ s| (* 3?) + 11 (If) + u If (4) 

In order to find the quantities (ST/c)y) and T_ shown in 

w 0 

equation (l) , two independent solutions of the boundary-layer 
equations are necessary. The first solution is for (c)T/dy) w 

with the assumption that T is independent of x and, with the 

following boundary conditions: 
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When 


J = 

u = 0 = v 

and 


and when 


T = T = constant 
w 


and 


7 = 6 , 


u = U 



The second solution is for T q with the assumption that 
(cHl/c3y) = 0 and T = T , which is an unknown quantity to he deter- 

"W* Q W 

mined hy solving the houndary -layer equations . 


Dimensionless Parameters 

The solution of equations (l) to (4) will involve the quantities 
u, v, x, y, p, \i, c , k, T, p, H, V, 0), cp, and 0. The 

quantities cp , and 0 are exponents of the temperature relations 
of viscosity, thermal conductivity, and specific heat, respectively. 
When the equations are put into dimensionless form, and solved, the 
previously, mentioned physical quantities appear as a group of dimen- 
sionless parameters, so that the solution can he expressed in the form 

'lKWu, Re, Pr, M, Eu, T^, A, to, cp, 0) = 0 (5) 


« 


4 
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where 


Ku = — = ISPusselt number 

K 


Re = — 2— = Reynolds number 


°-n w ^w 

Pr = ■ = Praruitl number 

Hr 


M = — = Mach number 
a 


Eu = g = Euler number 


T = -2. = temperature ratio 
w 


T e - T 6 
tt/zbJo 


= recovery factor 


= 

Ic 
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Hie Hub selt mmiber EV/V^ is a well-known heat -transfer param- 
eter that can he derived from equation (l) by rearranging and multi- 
plying "both sides of the equation hy the reference length x, thereby 
making equation (l) nondimensional. 

Hie Reynolds number TTp^/p^. is a well-known parameter for 
fluid -flow problems relating the inertia to the viscous forces. 

The Prandtl number c results from, making equation (4) 

rto md imftnfl l ortfl-1 r and depends entirely upon the properties of the fluid. 

The Mach number Tj/a and the temperature ratio Tg/T^ appear in 

the solution when the heat generated in the boundary layer is not 
negligible and when large temperature differences , producing changes 
in property values, exist between the fluid and the body. 

The Euler number is not as familar a grouping as the 

p s tj 2 /x 

preceding numbers. It measures the pressure gradient in a non d i m en- 
sional fashion. Although the symbol for Euler number has appeared as 
m in previous bound ary -layer investigations (references 5, .6, and 9), 
the symbol Eu is used in this report to conform, to the notation for 
other dimensionless quantities, such as Re, Pr, and Hu. When the 

velocity varies as x® 11 , the partial differentials in the numerator 
become total differentials and the Euler number reduces to the con- 
stant value Eu, as shown in references 5 and 6. 


The quantity T g is often related to Tg by the equation 

U 2 


T_ = T r + A _ 

© 6 2gJ c 


(7) 


where A is a dimensionless mmiber called the recovery factor. 


If large temperature differences exist between the gas and the 
solid surface, many of the physical quantities in the equations 
(c , fi., p, and k) vary by large amounts in the boundary layer. 

Jr 

In figure 2, the values of \x , k, and are plotted in a 

logarithmic scale against absolute temperature. The large varia- 
tions may not result in large changes in the Husselt number because 
seme of these variations may be compensatory. These curves can be 
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closely approximated "by straight lines on a logarithmic plot for mod' 
©rate ranges of temperature, so that it is possible to write 

v-k = (T/r f 

w V 

*/\. = (T/D w )'P 

y° PJ v " 

where CO, cp, and 3 are exponent a that will appear in equations (2) 
to (4) when the preceding substitutions are made for n , k, and 
c^. The usual substitution for p is made; that is, p varies 

inversely as the temperature and directly as the pressure. In the 
solution, the values of CO, cp, and 3 have an influence on the 
results when T r is markedly different from. 1. 



Effect of Variation of Dimensionless Parameters on 
Bus selt Number for Laminar Boundary Layer 
When solved for Bu, equation (5) takes the form 

Bu = ^ (Ee, Pr, M, Eu, T r , CO, cp, 3) (9) 

and 

A (Pr, M, Eu, co, cp, 3) (10) 

The evaluation of the functions and \|/^ presents a formi- 

dable mathematical problem. The BACA Lewis laboratory is now working 
out a numerical solution for the general case, where the gas proper- 
ties in the boundary layer vary with the temperature and large pres- 
sure gradients are present in the free stream at the same time. 
Calculations have already been made by several investigators and in 
this paper for some simplified cases by which the first-order effects 
of the parameters, acting independently, can be determined. The 
effects on the Buss el t number of the other nondlmensional numbers in 
equation (9) follow. 
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Effect of Reynolds number Re. - All the references dealing 'with 
heat transfer through a laminar boundary layer (references 1 to 9 and 
15); no matter ■which of the other parameters are included or emitted; 
agree that Hu . should be proportional to the square root of the 
Reynolds number (for example; reference 15 , pp . 626-636). Thus; 

Hu CD /s/Se” (ll) 

Effect of Prar^tl number Pr. - A simple and exact relation for 
the effect of the Prandtl number on the Husselt number is impossible 
because of the -way Pr enters into equation (4) when the equation 
is ttw/Ia 7-iortrj~?niAriA -T nn al - The effect has been calculated (reference 15; 
pp. 624 ana 626) in the case ■where M ~ 0, T r £ 1; and variations 

in c are not considered^ that is , p = 0 . It has been found that 

Jr 

with Eu * 0 for 0.7 < Pr <, 15 , the approximate relation is 

l/3 

Hu CD (Pr) ' (12) 

Similarly; with Eu = 1 for 0.6 < Pr < 2 (reference 15; p. 632); 

Hu CD (Pr)^* 4 (13) 

A similar relation for Eu = 2 presumably could be found, but as yet 
has not been determined. 

Effect of Euler number Eu. - The equations that give the effect 

of the Euler number on the Husselt number for the case when T ® 1. 

r 7 

14 «0, *r»a p =0; are given in references 5 and 6. Calculations 
are made in reference 5 for a series of Prandtl numbers from 0.6 to 
1.1 for Eu = 0 and Eu = 1 and for a series of values of Eu for 
Pr = 0.7 in reference 6. By utilizing the laminar velocity distri- 
butions published in reference 17 , Hu/^/E© may be calculated from 

the equation given in reference 6, which is 
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The results of the calculations are shown in figure 3, where 
-0.09 < Eu < 1.5 and 0.6 < Pr < 1.1. The lower value of Eu, where 
8u/dy is zero at the wall, gives the largest adverse pressure gra- 
dient theoretically possible in a laminar boundary layer with s 1. 

If this gradient is exceeded, reverse flow will begin ab the wal l and 
transition or separation will occur. When the value of Nu/V Re 
plotted in figure 3 is defined as E n , then 

F i >aai =w^ (15) 

where E.. is a function of both Eu and Pr. 
lam 

Effect of temperature ratio T r . - In the following paragraphs, 

the effect of T^ will be evaluated by keeping Ct>, op, and {3 of 

equation (8) constant, followed by the evaluation of Nusselt number 
changes due to variations in d>, cp , and £ . 

By utilizing figure 4, which is obtained from the results of 
reference 9, it is seen that at least for Eu = 0 the Uueselt number 
and the temperature ratio may be expressed in the form 

Hu OG (T r ) n (16) 

Then n can be evaluated from published results for the case when 
6=0 (constant specific heat) and Eu = 0 (no pressure gradient). 
The results of the evaluations where co and cp were assumed equal 
are shown in the following table: 


T r range 

Mach 

number 

Prandtl 

number 

0)= Cp 

n 

Eeference 

1 £ T r 4 

■ 0 

1 

0.76 

-0.089 

2 

1 < T <4 

0 

0.725 

.75 

-.073 

4 

1 < t£ < 4 

.9 

.7 

.7 

-.091 

3 

1 < T < 4 

5 

1 

.8 

-.087 

3 

1 ^.T <L 2 
r 

5 

1 

.8 

-.084 

3 
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CM 

a 


The values in the table indicate that neither changes in the Mach 
number nor in the Prandtl number have much effect on n in the ranges 
shown. It will be shown later that this Independence is not true of 
oj. The influence of to is not shown in the table because the values 
used in the references cover only a very small range. In practice , 
co cannot change much unless extreme temperature ranges are used or 
unless the fluid is changed. 

The wall temperature has been used throughout this investigation 
in evaluating ITu, Re , and Pr, all of which involve some of the 
fluid properties. If the free -stream temperature is used instead of 
the wall temperature for evaluating the fluid properties, the value 
of n will be quite different. This effect is shown in figure 4. 

The question immediately arises as to whether there exists some 
intermediate temperature that, if used for fluid-property evaluation, 
would cause n to become zero. 

It is suggested in reference 9 that same results of Crocco can 
be closely approximated by evaluating the fluid properties at a 
temperature T when 



(17) 


0 <M < 5 

Because in this correlation constant specific heat (P =0), zero 
pressure gradient (Eu = 0) , 0) = 0.75 = cp, and Pr = 0.725 are 

assumed, it is not recommended for other cases when these conditions 
are not fulfilled. 

The temperature ratio has another effect not directly related to 
equation (9) , the effect on the stability of the laminar boundary 
layer. A theoretical analysis is presented in reference 18 and the 
results are shown in figure 5, where the critical Reynolds number is 
plotted against the temperature ratio for Euler number of zero. The 
critical Reynolds number, which indicates the stability of the laminar 
boundary layer or its resistance to transition, is seen to increase 
for increasing T^. . It is experimentally shown in reference 19 that 

the Reynolds number Of transition decreases when T becomes less 
than 1. r 

Effect of viscosity-temperature exponent to. - In going from one 
gas to another, the value of CO may change considerably. It is shown 
in reference 20 that co ranges from 0.647 for helium at low tempera- 
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tures to 0.98 for carbon dioxide. Tlius a range of Tallies of CO could 
conceivably be useful. Reference 4. contains such values with M =lj 
Pr = 0.725, gp= c p , P = 0 , and Eu = 0 (presented herein as fig • 6 ) , 
wherein c f g, a quantity proportional to Busselt number divided 

is plotted against CO. This figure shows that when 

= 4, a 


L,6 


i>y 

T r -1, 

change of 
decreases IJu 


Bu^ qA/ _ ^®l~ 5 is a1moa ^ i^depandejn-t of CO. If T^ 

(0 from. 0.7 t o 0.8 (plausible range when the fluid is air) 


L,6 


A / 


Ee 


L,5 


frcffli 0.65 to 1.0 decreases 
ence value of 




L,5 


about 3 percent , whereas a change in co 
■t. g/^Re^ ~ about 8 percent , the refer- 
being at T = 1. 


Effect of thermal conductivity - temperature exponent co* - The 
separate effect of cp. has not been published insofar as is known. In 
the references it is assumed to be the same as £0. 


Effect of specific heat - temperature exponent P . - The effect 
of the exponent p has been computed using the results of refer- 
ence 4 and is shown in figure 7 with K, a quantity proportional to 
Efu^ g/^Re^, ™ , plotted against p. This figure shows that for 

T_ = 1, HiL fe /./Re T is independent of P . The figure also i llus- 


t - - 

trates the fact that 


Ku l, g/^Re^ 5 decreases for increasing T^. 

At P =0.2, a plausible value fo r air, K is reduced 10 p ercent for 
T r = 4.0, ^.inasmuch aa MK, 

reduced 10 percent from, the value at T =1.0. 


is also 


It is therefore apparent that the large changes in the physical 
properties of the fluid through the boundary layer due to large tem- 
perature variations produced several effects. Some of these effects 
tended to reduce Ru//s/Re, whereas other effects tended to increase 
Bu/a/R©. The net result of all these changes has not been completely 
calculated, but In view of the compensations noted, the net effect 
may be much less than the changes in the individual properties would 
indicate . 


Effect of Mach number. - Calculations are presented in refer- 
ences 3 and 4 in which a range of Mach numbers was used. From an 
analysis of the curves based on these calculations. It Is concluded 
that for Mach numbers not exceeding 4 an approximate relation is 



1346 


BACA m E5QF02 


13 


Tlie coefficient C in agnation (18) depends on the other param- 
eters of equation (9) , hut does not undergo large changes in the range 
of values usually encountered. The calculated variation caused by 
changing Mach number frcaa 0 to 2 is shown in figure 4. When 

(0 = 0.75 = qp , Pr = 0.725, T^ = 2, and {3 = 0 = Eu, a calculation 

utilizing the results of reference 4 yields the value 0.00374 for C. 
This value was altered very little hy changing from. 1 to 4. 

The curves of reference 3 indicate a value of 0.00341 for C, where 
T r = 1, Pr = 1, a) = 0.8 = cp, and Eu = 0 = p. Both of these num- 
bers are so small that for subsonic flow, the Mach number effect is 

quite inappreciable, at least for Eu = 0. When M = 4, the effect 
on Hu is to reduce it about 6 percent from the value at M - 0. 

Correlation equation for laminar flow. - For 1 1 T < 4 and 

0 M <.2, the combined effect on the Husselt number for air as the 
gas of T r , M, co, and P is of the order of a 5- or 10-percent 

change from the value at T =■ 1, M = 0 (for air, P =0.2 and 

CO = 0.7). Xn view of the compensations on Hu/a/B© of the different 
factors involved, an approximation has been made for purposes of cor- 
relation for the heat transfer frcan a laminar boundary layer. For 
this approximation, equation (9) takes the form 

Hu =^ 3 , (Be, Pr, Eu) (18a) 

Figure 3 was put in the form of figure 8 by using the results of 
figure 3 for Prandtl number of 0.7 and results of reference 6 for 

Eu > 1.5. For the range 0 <l Eu < 2, the results shown in figure 8 

check figure 3 to less than. 3 percent. Thus the following equation 
replaces equation (15) far 0.6 < Pr < 1 and 0 < Eu< 2: 

Nu/(Pr) l/3 = ijJS (19) 


where 


=*W ^ 1/3 

Equation (19 ) will be the correlation equation for laminar flow with 

F given in figure 8. 
lam 
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Equations of Turbulent Boundary Layer 

When the boundary layer becomes turbulent , the differential 
equations (2) to (4) can no longer be applied because the motion is 
not steady and the mean motion does not satisfy the equation of 
Stokes. In spite of much simplification, neither the dynamic nor 
the ther mal problems are at present amenable to solution in their 
more exact form. The flow is often regarded as consisting of a 
mean flow and a superposed fluctuating motion. By use of this 
concept, equations similar to equations (2) and (4) have been pre- 
sented wherein the velocities and the temperatures are time averages 
(reference 8) that neglect pressure -gradient and dissipation terms 
as follows: 

u S +v 57- + 5?(p§? +e M^I 

t h St . - dr\ d ( ul 8t 8t\ 
s ( u S + T &) - 5? + C H 557 

where is the so-called eddy diffusivlty for momentum and 

is the so-called eddy diffusivlty for heat. Thus formally, the 
laminar solutions for no pressure gradient and, no dissipation could 
be used by substituting ^ + p/p for p/p and € E + p/p Pr for 

p/p Pr. This solution cannot be carried out so simply, however, 
because the values of the eddy diffusivities vary from place to place. 
Consequently, equations (20) and (2l) have been used to derive Reynolds 
analogy. This derivation is carried out by postulating that 
dp/8cx: =0, Pr = 1, and Equations (20) and (2l) then 

become identical in form, the_only difference being the substitution 
of T in equation (2l) for u in equation (20) as the dependent 
variable. It can then be shown (reference 8) that 


( 20 ) 

( 21 ) 


H 

c p ,w p w U 



( 22 ) 


where f is the friction coefficient equal to — =» — 

P q V 2 /2 

The postulate that Pr = 1 holds approximately for all gases. 

A number of formulas have been proposed for correcting equation (22) 
for cases where Pr / 1. The simplest formula, given in reference 21, 
is 
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H 




(Hr) 


2/3 


_ f 


(23) 


It is shown in reference 8 that equations containing more accurate 
ftnri more complex correction factors do not differ markedly from, 
equation ( 23") in the r ang e 0.5 < Pr < 10. "When equation (23) is 
solved for Hit, the result is 


Nu = | He (Pr) 


1/3 


(24) 


The correlation for Hu can then he completed hy substituting 
an f correlation in equation ( 24) . One of the" simplest of these 
correlations Is given In reference 8 and Is 

£ = 0.0296 (Re)’ 0 * 2 (25) 

2 

Substitution of this value in equation (24) yields 

Hu = 0.0296 (Re) 0 * 8 (Pr) 1 ^ 3 (26) 

Reference 9 presents equation (26), compares It with two more 
elaborate correlations, and concludes that equation (26) Is adequate 
for .gases , though possibly not for liquids with extreme values of 
Pr. Accordingly, equation (26) will be used as the correlation 
equation for turbulent boundary layers even though the pres sure - 
gradient effect has not been included. 

In the laminar boundary layer. It was found t h at the wall tem- 
perature was suitable for evaluating the gas properties p, c^, p, 

and k. Although the evaluation temperature for a turbulent boundary 
layer cannot be obtained fTcsn. this analysis, a temperature equal to 
the surface temperature is suggested, considering that the laminar 
sublayer in the turbulent boundary layer may be controlling the heat 
transfer. Also, good correlation for turbulent flow In pipes Is 
obtained in reference 22 when the gas properties are evaluated at the 
wall temperature. The wall temperature will therefore be used for 
evaluating the gas properties In both laminar and turbulent boundary 
layers . 
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Heat Transfer with Both Laminar and Turbulent Boundary layers 

As the solid surface increases in length, the Reynolds number 
■becomes -larger and the laminar layer ultimately becomes unstable, 
changing to a turbulent boundary layer. The problem then arises of 
computing an average Uuseelt number for the whole body, from the 
stagnation point where x = 0 to the trailing edge where x - L. 

If the transition point is at |L, where the l ami nar boundary layer 
changes to a turbulent boundary layer, then 

(27) 

transfer is to integrate H 
layers , 


dx = H 1 (28) 


The heat transfer over the whole surface will be obtained when 
both laminar and turbulent boundary layers are considered to exist 
on the surface. 

Heat transfer through laminar boundary layer.' - The correlation 
equation for the laminar boundary layer is given as equation (19) , 
which may be written as 



The method for obtaining the total heat 
for both laminar and turbulent boundary 



H =r- -1/2 

^173- w\| 


(29) 


In order to obtain the heat transferred, an integration over the 

laminar layer is necessary. The assumption of constant p will be 

v 

made, which is reasonable for low Mach numbers. In addition, the wall 
temperature and 3i *i will be assumed constant over the laminar 

portion. , Then integration of equation (29) yields 
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Xt is stated in references 5 and 6 that if Eu is to he con' 
stant, the relation "between U and x is 


U 


°i x 


Eu 


(31) 


where Eu is the Euler number. Substituting equation (3l) in 
equation (30) and integrating yield 


k v (Er) 


173 


rt L 


E dx 


2 E Ip c", 
lam / w 1 

Eu+1 ^ H T 


Eu+1 


2 

«L) 


(32) 


Heat transfer through turbulent boundary layer. - Xhe corre- 
lation equation for the turbulent boundary layer is given as 
equation (26), which may be written as 


i^(Pr) 1/5 



(33) 


Proceeding in a manner similar to equations (29) to (32), with the 
integration now from |L to L, equation (33) becomes (neglecting 
pressure -gradient effects on the local heat-transfer coefficient) 



Summation of heat transfer over surface length L. - "Where both 
laminar and turbulent boundary layers exist over tile surface, use of 
equations (28), (32), a n d (34) w i ll yield an average heat -transfer 
coefficient, thus 
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HL 

y^j 175 


2F n /Z3T5T 
l«rn J M W 1 

Eu+1^/ M w 


Em-1 

2 

ftU + 


0.037 [ C l p v\ °* 8 l 0-8(Eu+D f .0.8(Eu+lj| 

Eu+1 U / *- -J 


(35) 




l/3 


The left side of equation (35) is lTu^/(Pr) , the ordinate custom- 


arily used in experimental correlations. The right side is usually 
given as a function of a Reynolds number. A suitable Reynolds number 
for use vith UL is 


Up L 

V 


(36) 


where U is the average velocity over the solid surface. It is 
necessary to express c-^ in terns of TJ. 

The defining equation for TJ is 


UL = I TJ dx 

Jo 


(37) 


Eu 


Substitution in equation ( 37 ) of U - c^x and integration yields 


o (38) 

1 Eu ' ' 

Xi * 

Substitution of this value of c^ in equation (35) yields 


Eu 


(Pr) 


I7i 


^laa r=T t 2 
V®ST^ Ee L 5 


Eu+1 


0.037 — 0.8 

^2 Be L 


(Eu+1) 


1-1 


0.8(Eu+l) 


(39) 


The reference length used in the Kusselt and Reynolds numbers is 
sometimes different from, the L used here. For turbine -blade shapes f 
a common reference length is perimeter /it } which reduces to the 
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diameter for a circular cylinder. In equation (39), perimeter/* 
will be used by substituting 2L/* for I in Hu^ and. in Be^« 

The coefficient 2 In the first tern then becomes 2 \/2/jt or 
1.596, -whereas the coefficient 0.037 in the second term becomes 

0.037 (2/*) 0,2 or 0.0338. 


l/3 

If values are chosen for Eu and £, then Nu /(Er) 7 Is a 

E 

function of Ee^. Even when plotted on a logarithmic scale, the 


function will still be slightly curved because of the 0.5 and 0.8 
powers of the Beynolds number in equation (39). Because heat "transfer 
correlations are usually presented by drawing a straight line through 
the observed points plotted on a logarithmic scale , for purposes of 
comparison equation (39) will be put into the approximate form 


Hu /( Er ) l/3 = F (He f 

p p 


(40) 


with perimeter/ jr being the reference length. 

Equation (39) can be put into the form of equation (40) by 
malrTng the exact equation (39) coincide with the approximate 
equation (40) at two points near the ends of the observed range of 
experimental data. In this case, the two points chosen will be 
- 5 

E©p = 2 x 10 and 2 x 10 • The greatest difference between the 

exact and approximate curves will then be about 2 percent or less 
4- 5 

for 10 < Be < 4 x 10 , increasing to about 6 percent when Ee 

^ 6 p 

is extended to 10 . 

Values of z and F are shown in figure 9. They are plotted 
as ordinates against the transition ratio £ as abscissa. Four 
values of Eu are shown, ranging from -0.09 to 2. The effect of 
Eu in figure 9 is less than that shown in figure 8 because Eu 
appears in the denominator of the first term, of equation (39) and 

thus partly cancels the effect of F„ 

lam 


By using figure 9 and equation (40) , curves were determined for 
*3i^f®^sut ^values of £ and Eu, as shown in figure 10, where 

(^ ) is plotted against Be^. As would be expected, the 

lowest heat transfer results frcaa a completely laminar boundary 
layer, £ = 1; whereas the highest heat transfer is obtained for an 
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entire turbulent boundary layer, | = 0. In practice, some of the 
curves would not extend into all the Reynolds number ranges shown, 
but the extension has been made in figure 10 to illustrate the 
theoretical trends. It is seen from the slope of these theoretical 
curves why _ experimental results, when plotted in this fashion, often- 
times yield slopes different from 0.5 or 0.8, the theoretical values 
expected with either completely laminar or turbulent boundary layers, 
respectively. 


APPLICATION OP THE GET TO EXPERIMENTAL RESULTS 

In the following section, the theory will be compared with exper- 
imental results obtained from cylinders, an airfoil, and turbine -blade 
cascades. On these surfaces, Eu varies locally but average values 
will be used except for stagnation points and the airfoil where local 
heat -transfer coefficients have been measured and can be compared with 
theory. 


Cylinders 

Equation (19), when written in the form 

■\ fx /p IbsT 

? (41) 

w v 

can be used to calculate a Efusselt number for the stagnation point of 
a cylinder. Near this point, it has been shown (reference 23) that 


U = 3.63 ~~ (42) 

where D is the cylinder diameter and IT the upstream velocity 

, of the fluid. Substitution of equation (42) in equation (41), 
cancellation of the x*s, and multiplication of both sides by D 
yield 


m 

K 


1.905 P. 


lam 


(Pr) 


l/3 


P UJD 
w 1 


w 


(43) 


Por this case_ 2 _ Eu = 1 (reference 15, p. 631) and utilizing fig- 
ure ( 8 ) for P_ with Pr =0.7, equation ( 43) becomes 

JmQHH 



HACA EM E50F02 


21 


H 

.Ctf 

ff»- 

o> 



or 



(44) 


(45) 


Equation (44) is nearly the same as that given in reference 15 
(p. 632) . The difference in the numerical coefficient (l.Ol) in the 
reference is due to the use of a higher Frandtl number (Pr = 0.733) 
and a theoretical velocity distribution (U = 4 U^x/d) . 

Experimental measurements of the Husselt number at the stagnation 
point of a cylinder are given in reference 24 and are plotted in fig- 
ure 11 for three cylinders of diameters 1.27, 2.54, and 6.35 inches. 
The theoretical curve is shown also and agrees very well with the 
data. 


The average Husselt numbers over the cylinders are shown in fig- 
ure 12, together with the theoretical curve from equation (40). The 
Reynolds number Ee^ containing the upstream velocity has been 

converted to the Reynolds number Re^., which contains U, the 


average velocity around the cylinder. The average velocity and the 
dimensionless number Eu were found from, the 'pressure distributions 
given in reference 24 through use of Bernoulli 1 s equation and the 
definition of the Euler number. The value of £ was obtained by 
noting the point of minimum Husselt number in the plot of Hu 
against the angle measured from the stagnation point. Thus, £ was 
found to be 82°/l80° or 0.455. The agreement between the experiment 
and the theory is surprisingly good in view of the fact that separa- 
tion has occurred on the cylinder and no account of separation was 
taken in the theoretical development. 


Airfoil 

Reference 25 gives local heat -transfer coefficients for clear - 
air conditions obtained in flight investigations of a HACA 65,2-016 
symmetrical airfoil with an 8 -foot chord. An equivalent diameter 
at the leading edge was obtained herein using the coordinate sys- 
tem for the airfoil given in reference 26. From this equiv- 
alent diameter and the data of reference 25, it was possible to 

1/3 

determine HUp/(Pr) and Re^ at the stagnation point. Because 
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only one flight speed -was used for the clear-air condition, only one 
point was obtained, which is shown in figure 11 and agrees well with 
the theoretical curve and experimental results from cylinders. 

Using the experimental results of reference 25, it was also 
possible to determine localjiusselt , Reynolds, Prandtl, and Euler 
numbers. The coefficient 3T could then be determined loc all y 

from the Euler number and a theoretical curve drawn fear comparison 
with the experimental results. The comparison is made in figure 13, 
equation (19) being for the laminar case and equation (26) for the 
turbulent case. The curve for a flat plate is given for the laminar 
range to show the increase resulting from inclusion of Euler number 
effects in equation (19). A scale is also given for the variation 
of x/L with Re. It is seen that equation (19) agrees well with 
the data in the laminar range until very close to the leading edge. 
Equation (26) falls somewhat below the experimental results. 

„ A/3 

An average value of Nu /(Pr) was calculated using equa- 

p 

tion (40) and figure 9. The average Euler number was obtained from 
the pressure distribution and the transition ratio was obtained from 
the plot of local heat -transfer coefficient against x/L. The tran- 
sition point was taken as the value of x/L where the heat -transfer 
coefficient shows a sharp increase. By using this method, the calcu- 

l/3 

culated value of Hu /(Pr) ' was about 9 percent less than the 

P 

experimental value. 


CO 


a 


rH 


Turbine -Blade Cascades 


In order to determine average outside coefficients for turbine 
blades, various sets of heat -transfer measurements have been as- 
sembled. The pertinent facts about these experiments are shown in 
table I. It is evident from the table that various pressures, tem- 
peratures, velocities, and lengths were used by the different authors 
in their correlations. In order to put these results on a uniform 
and comparable basis, the gas properties were all reduced to those 
at blade temperatures , the velocities were reduced to average values 
around the blade periphery outside the boundary layer, and the refer- 
ence length was reduced to perimeter /« . The velocity *reduction was 
accomplished by using the state of the gas at the inlet to the blade 
cascade, calculating the velocity distribution around the blades 
through use of the stream-filament method, and then utilizing the 
experimental curves and the data where available to obtain 
l/3 

HUp/(Pr) . The details of these calculations are given In refer- 
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ence 14. Scan© of tile experimental correlations resulting are shown 
in figure 14. 

The comparisons "between the correlation curves given in figure 14 
and. the theory are made in figures 15 to 17. The method used to 
obtain the theoretical curves is given in appendix B. The comparisons 
indicate better agreement of the theory and the exper ime nt when f is 
determined by the condition that Eu= 0 rather than by the condition 
that Eu = -0.09. An exception is. the case of reference 13 (fig* 17) , 
where better agreement is obtained for the Eu = -0.09 condition. 

The deviation between the f determinations in this example , however, 
is rather small. 


SUMMAET OB' KSSUITS 

From a theoretical analysis of heat transfer from, the boundary- 
layer aspect, the following results were obtained: 

1. Other dimensionless numbers in addition to the Reynolds and 

Erandtl numbers are needed for an accurate evaluation of the Busselt 
number. Two of the most influential numbers are: the Euler number, 

a measure of the pressure gradient; and the transition ratio, which 
measures the amount of laminar boundary layer present on a surface. 

2. The Mach number appeared to have a smaller effect, at least 
for Mach numbers less than 2 and temperature ratios between 1 and 4, 
when the fluid was air. The viscosity- and specific heat' - tem- 
perature exponents showed varying influence , depending on the tem- 
perature ratio. 

3. An equation for the average heat transfer was derived for 
constant wall temperature to Include the Euler number and transition- 
ratio effects. 

4. Theoretical curves based on this analysis agreed with experi- 
mental results obtained from cylinders, a symmetrical airfoil, and 
turbine -blade cascades. 


Lewis Flight Propulsion Laboratory, 

National Advisory Committee for Aeronautics, 
Cleveland, Ohio. 
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SYMBOLS 

Tlie following symbols are used, in this report: 
speed, of sound, in gas 

factor in correlation equation (18), Eu co (l - C M ) 
drag coefficient due to akin friction 

specif io heat at constant pressure at point x,y in boundary- 
layer 

Eu 

factor in U = c^x 
diameter of cylinder 

2L 

characteristic length, — 

Euler number, - 

P 5 xr/x 


F 


‘lam 


E 


St 


mean coefficient , 

(Pr) 7 


coefficient for laminar flow, 


Eu 


F. 


lam 


mean coefficient. 


(Br) 


la m 

173 


friction coefficient, 


v 


P*U 2 /2 


s 


nondime nsional stream function 




Eu+1 
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S 

H 

H 

h 

J 


fc 

L 

L tr 

M 

Hu 


acceleration due to gravity 
outside heat “transfer coefficient 
mean outside heat -transfer coefficient 
enthalpy 

mechanical equivalent of heat 

h -iu 

frcm reference 4, K = 7 = — ■ [— 

(T -T-Jc 7 
v v 6 p,S 

thermal conductivity at point x, y in boundary layer 

total length of surface, measured along surface from stagnation 
point to trailing edge 

value of x at transition point 
Mach number, u/a 

local Husselt number based on x, Hx/k^ 
local Husselt number based on D, HD/k 


mean Husselt number based on L, HL/k 


■w 


Hu mean Husselt number based on ^ HD /k 

P it o v 

n exponent of Husselt - temperature -ratio relation, Hu CD (T^) 

Er Erandtl number, ^ 

p static pressure in boundary layer, or at its edge because 

3p/^y is assumed zero 

q stream function of flov in boundary layer 

E universal gas constant 

He local Eeynolds number based on x. Up x/x 

V V 
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E© critical Reynolds number below which all oscillations in 
cr laminar ‘boundary layer are damped 

Re_ Reynolds number "based on D, XL o D /u. 

D 7 1 v v 

Re_ mean Reynolds number based on L, Up L/p 
L w w 


Re 

T 

T 

T 

P 


T 

] 

U 

U 



u 

u 

V 

V 
X 


J 

Z 



distance normal to surface 

exponent of Reynolds number (equation (40)) 


1346 



9^21 


NACA EM E50EFQ2 


27 


0 

5 

C H 

£ M 

A 

V- 

V 

i 

p 

T 

V 

cp 


* 

*a 

L 

>1/ 


CO 


exponent in specific heat - temperature relation, 

O Jo = (T/DJ P 

P p,v V 

outer edge of boundary layer 


eddy diffusiYity for heat 
eddy diffusiYity for momentum. 


nondimensional boundary -layer coordinate, 

T - 2L 


recovery factor , 


TJ 2 / 2 


absolute Yiscosity of gas at point x, y 
kinematic Yiscosity (u/p) 


transition ratio. 


tr 



in boundary layer 


density of gas at point x,y in boundary layer 
shear stress, (jiSu/^y)^ 

exponent in thermal conductivity -temperature relation, 

k/k = (tA 

v -w 

function, equation (5) 
function, equation (9) 
function, equation (10) 
function, equation (l8a) 

exponent in Yiscosity-temperature relation, \i/\i = (t/t )° ) 

v v 
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Subscripts: 

D based on diameter 

e effective 

L bdsed on total length 

lam laminar 

P pressure surface 

p based on perimeter 

S suction surface 

SP suction and pressure surface 
w vail condition 

8 main stream condition 
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APPENDIX B 


METHOD 05* EVALOATIFG- TRANSITION - RATIO AHD AYERACT EjJLER HUMBER 

In order to obtain the theoretical curves for comparison with 
experimental correlations shown in figure 14, the chor&wise velocity 
distributions given in reference 14 were used to determine the point 
of transition and the average value of the Euler number, inasmuch as 

Eu = - (HL) 

PJJ /X 


and from Bernoulli’s equation 


Equation (Bl) becomes 


dp TT dU 
- = pTJ -rrr 


EU 

TT dx 


By using finite differences, dU/dx was calculated from the 
aforementioned velocity distribution and combined with the other 
factors in equation (B2) to obtain Eu« The resultant plot of Eu 
against x/L is shown in figure 18 using the velocity distribution 
of reference 12 given in reference 14. As shown in figure 18, the 
transition point was obtained at x/L, where Eu -0 or Eu = -0.09 
for both the suction and pressure surfaces. The value Eu = -0.09 
represents the condition where the adverse pressure gradient is so 
large that reverse flow at the wall is imminent. The value Eu = 0 
represents the limit of the adverse pressure region. At any value 
of Eu less than Eu =0, the boundary-layer velocity profile has 
a point of inflection, indicating an unstable velocity regime. 

A weighted average of the transition ratio was determined using 
the lengths of the suction and pressure surfaces. 


+ 

+ I * 
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The Euler number -was determined "by obtaining the integrated 
average value of Eu to the point of transition for both the 
suction and pressure surfaces. A weighted average was also used 
for the Euler number. 



L Eu + L Eu 
c S S P P 


+ L 
P 


(B4) 


Using the values of % and Eu^ calculated by equations (B3) 

and (B4) , figure 8 was utilized to determine Z and' F and these 
values substituted in equation (40) to give the theoretical curve. 
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TABLE I - SU&MARI OF VARIABLES USED IN BEAT-TRANSFER EXPERIMENTS ON OOTSH2E OF TURBINE BLADES 


Investigator 
Reference number 

pollmann 

11 

Meyer 

12 

Andrews, Bradley 
13 

Temperature ratio, T r 

0.8 

0.9 

2-1 

Evaluation temperature for: 
vlsoosity 

Inlet stream 

Film 1 

Blade 

density 

Inlet stream 

Inlet stream 

Mean 2 

thermal conductivity 

Inlet stream 

Film . 1 

Blade 

Pressure for density 

Inlet 

Inlet 

Mean 2 

Velocity 

Inlet 

Inlet 

Exit 

Length 

Perimeter/n 

Perimeter/if 

Chord 

Chord, in. 

3.94 

0.680 

1.0 

Number of blades 

3 

6 

5 

Solidity 3 

1.47 

1.92 

1.61 

Method of heating or cooling 

Electrically 
heated blade 

Electrically 
heated blades 

Hot gas stream with 
water-cooled blade 

Remarks 

Only middle 
blade heated 

Blades heated 
by conduction 

Only middle blade 
cooled 


Film temperature = l/2 (stream static temperature + wall temperature ) . 
^Mean density « pA^^an 

p •= l /2 (inlet pressure + exit pressure) from cascade 

T = 1/2 (T + T ) , T - T S + 0.85 U 2 / 2 gJc 
■7 mean ' e w e * 

■ (Solidity - chord /spacing 
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(a) Fluid properties based on. wall temperature. 
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Temperature ratio, T r 

(b ) Fluid properties baaed on stream temperature. 

Figure 4. — Theoretical effect on heat transfer of temperature ratio and Mach 
number for gaa properties at wall or stream temperature (reference 9). 

Pr = 0.725; to — 0.75 = cp; p = 0; Eu = 0. 



Temperature ratio, T r 

Figure 5* — Theoretical effect of temperature ratio 
on stability of laminar boundary layer (reference 
18). M = 0.7; co = 0.76; Eu = 0; Pr = 1.0. 




1 . 0 , 



Figure 6. — Theoretical effect of co on heat transfer 
with laminar boundary layer (reference 4 0« M = 1.0; 
Pr * 0.725; <-0 = (p; ft = 0; E u = 0; }X Oz T^; 

Nu L,s// Re L,5 ocCf/Re^. 
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Transition ratio, £ - L$ v /L 


Figure 9. - Effect of transition and 
correlation equation Hu^/CPr) 1 / 3 


Euler number on 
= F(Re p )Z. 


0 













Diameter 

(in.) 

6 . 35 ^ 

2,54 ■ reference 24 
1 , 27 , 

- Theoretical curve, 
equation (40) 




Figure 12, - Comparison of experimental data for average Nusselt number from 
cylinders with theoretical curve for Eu “0.55 and f - 0.455 (data from 
reference 24), 
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Figure 13. - Comparison of local heat-transfer obtained theoretically and from experiment of 

reference 25. 













Figure 15* - Comparison of theory with experimental data of reference 11 

Impulse blade j temperature ratio, 0*8. 




Figure 16# - Comparison of theory with experimental data of reference 12 

Impulse blade; temperature ratio* 0,9, 




Figure 17. - Comparison of theory with experimental data of reference 15 

Reaction blade; temperature ratio. 1,22. 
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Suction surface x/L Pressure surface 

Figure 18. — Variation of local Euler number for reference 12. 
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